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ABSTRACT 

A  natural  circulation  loop  wi.th  water  as  the 
circulating  fluid  was  studied  for  a  range  of  operation 
covering  two-phase  flow.  Periodic  oscillations  of  the 
flow  rate  and  fluid  temperature  occur  even  with  constant 
heat input  and  constant  cooling  water  properties  for  the 
heat  exchangers.  Several  conclusions  concerning  the 
stability  of  operation  are  given.  Use  is  made  of  the 
theoretical  analysis  of  an  open-ended  system,  and  an 
analogue  computer.  For  use  in  a  more  detailed  numerical 
analysis,  the  equation  of  motion,  the  continuity  equation. 


and  the  energy  equation  are 
phase  flow  model. 


Dresented.. 


for  a  transient  two- 


The  Oscillatory' Behavior  of  a 

Two-Phase  Natural  Circulation  Loop 

by 

Eugene  H.  Wissler,  H.  S.  Isbin,  and  N.  R.  Amundson 

Introduction 

A  program  of  study  of  the  transient  operation  of 
natural  circulation  loops  has  been  underway  at  the  University 
of  Minnesota  (l),  and  this  paper  is  concerned  with  the 
oscillatory  behavior  of  a  two-phase  natural  circulation 
loop.  These  studies  are  of  interest  for  the  emergency 
cooling  of  nuclear  reactors  and  in  the  design  of  boiling 
water  reactors.  The  literature  survey  pertaining  to  the 
transient  operation  of  a  natural  circulation  loop  is 
given  by  Alstad,  Isbin,  Amundson  and  Silvers  (l),  and  a 
survey  on  two-phase  flow  is  given  by  Isbin,  Moen  and  Mosher  (2). 

Experimental  Loop 

Figure  1  is  a  schematic  diagram  of  the  natural 
convection  loop  which  was  studied.  The  loop  was  constructed 
primarily  of  l6  gage,  1  inch  O.D.  (0.872  inch  I.D.)  hard 
drawn  brass  tubing.  The  major  features  of  the  loop  are 
described  in  reference  (l).  During  a  natural  circulation 
run,  the  flowrators  were  by -passed  and  only  the  electro¬ 
magnetic  flowmeter  was  used.  Normally  the  surge  tank. 
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E^,  was  not  used  for  two-phase  natural  circulation  runs. 

For  those  runs  In  which  the  pressure  at  one  point  in  the 
loop  was  held  constant,  the  gate  valve  between  the  surge 
tank.  Eg,  and  the  loop  was  opened;  for  the  constant  volume 
run,  the  gate  valve  was  closed.  A  heater  and  pump  were 
installed  to  maintain  the  cooling  water  supply  at  5  gpm  and 


at  130 °F. 


General  Equations 


Theoretical  Analysis 

The  continuity  equation,  the  equation 


™  motion  and  the  energy  equation  for 

a  viscous  fluid  flowing  in  a  region  of  general  geometry 
were  formulated  for  a  phase  having  continuous  properties. 

A  similar  set  of  equations  were  derived  for  flow  across  a 
surface  of  discontinuity.  The  combination  of  these  two 
sets  of  equations  permits  one  to  write  the  equations  for 
the  two -phase  flow.  An  annular  flow  model  was  selected 
to  illustrate  some  of  the  essential  properties  of  two-phase 
flow.  The  model  is  sufficiently  simple  to  permit  attempting 
a  numerical  solution.  Each  phase  is  assumed  to  flow 
through  a  well  defined  cross-sectional  area  with  a  uniform 
velocity  with  the  reservation  that  the  liquid  velocity  at 
the  wall  must  equal  zero.  Further,  the  steam  and  water 
phases  are  assumed  to  be  in  equilibrium;  that  is,  p  and 
T  correspond  to  the  saturation  pressure  and  temperature. 
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For  these  conditions,  the  continuity  equation  may  be  written 
as 


*#--£  (1) 
The  equation  of  motion  for  upward  flow  through  a  pipe  of 
constant  diameter  becomes 


1_^W  1  3«1V) 

Sc  «  Sc 


(2) 


Finally,  the  energy  equation  is 


.  a  (oh)  ,  'd(wl|lgH)  a  Ha  _  0,  +S_  r(k  a  +kA)^i  (3) 

A  ot  +  Ox  -  J  -ot  -  Q  +  *^KWAw  +  KsAs'*’bxJ 

where  2vr  k  (2^)  =  Q*  =  rate  at  which  heat  is  con- 

p  w  or  r 

p  ducted  through  the  pipe. 


The  terms  and  $2  represent  the  ratio  of  the  true 
rate  of  transport  to  the  rate  of  transport  of  the  mean 
flow  for  momentum  and  energy  respectively.  If  both  phases 
have  the  same  linear  velocity,  0^  and  are  unity* 

The  three  functions  <J>2  and  Ff  were  deter¬ 

mined  experimentally  from  steady  state  data.  It  was  found 
that  and  <j>2  could  be  correlated  as  functions  of 

p  alone,  and  that 

Ff  =  afpJVJ1'79 


where  a(p)  is  a  function  of  p.  Figures  2,  3  and  h 
illustrate  the  variation  of  <J)g  and  a(p)  with  p. 


(4) 
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Equations  (1)  through  (4)  have  been  applied  to  a  natural 
circulation  loop  in  the  form  of  finite  difference  equations. 
The  discussion  of  these  equations  is  not  included  in  this 
paper  for  the  numerical  calculations  using  the  SEA.C 
(National  Bureau  of  Standard* 3  digital  computer.  Standard 
Eastern  Automatic  Computer)  have  not  been  successfully 
completed . 

An  insight  on  the  factors  which 

Stability  Analyses, 

determine  the  stability  of  a 

natural  convection  system  is  gained  through  the  analysis 

of  an  open-ended  loop,  such  as  shown  in  Figure  5.  The 

fluid  entering  the  heater  always  has  the  constant  temperature 

T  ,  and  the  velocity  of  the  stream  is  fixed  by  the  density 
oi 

difference  between  the  hot  and  cold  leg.  For  any  constant 
heat  input,  one  may  define  a  state  of  equilibrium  in  which 
the  difference  in  weight  of  the  two  legs  is  just  equal  to 
the  frictional  resistance  to  flow.  Under  certain  conditions 
this  system  may  be  unstable;  that  is  a  small  deviation 
from  the  equilibrium  temperature  distribution  or  the 
equilibrium  velocity  may  be  propagated  in  space  or  time 
with  increasing  amplitude.  For  example,  if  a  velocity 
perturbation  of  the  form,  Spinet ,  is  present,  under 
what  conditions  could  the  temperature  perturbation  which  is 
generated  cause  the  driving  force  that  sustains  the  velocity 
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dlsturbanee?  The  problem  ha3  been  treated  analytically 
for  a  one-phase  fluid  (3)»  and  the  following  conclusions 
were  obtained; 

a)  If  an  oscillatory  flow  rate  is  to  be  possible, 
the  force  cannot  be  generated  in  the  heater;  it  must  be 
generated  in  the  vertical  riser- 

b)  The  product  of  the  coefficient  of  expansion 
of  the  fluid  and  the  vertical  height  of  the  riser  must 
exceed  a  certain  value  (defined  by  an  analytical  expression) 
if  the  velocity  perturbation  is  to  be  sustained. 

c)  The  period  of  oscillation  will  be  approxi¬ 
mately  equal  to  the  residence  time  of  the  fluid  in  the 
heater  and  the  vertical  riser. 

In  order  to  predict  the  period  for  a  closed 
loop,  a  problem  was  solved  on  a  Reeves  Electronic 
Analogue  Computer.  The  number  of  non-linear  terms  was 
limited  to  16,  and  the  model  used  was  necessarily  a 
simple  one.  For  such  a  problem,  the  period  of  the  oscilla¬ 
tions  should  be  meaningful  even  though  the  wave  shape  Is 
not  correct,  if  the  period  is  a  function  of  only  the 
geometry  and  the  mean  velocity.  The  computer  could  handle 
only  ten  subdivisions,  with  the  driving  force  expressed  as 
a  linear  function  of  all  ten  enthalpies,  and  the  frictional 
force  as  a  quadratic  function  of  the  velocity.  The  number 
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of  available  summing  amplifiers  limited  the  problem  to 
the  case  in  which  boiling  occurs  only  at  the  top  sub¬ 
division  of  the  vertical  riser.  The  equations  solved 
were  of  the  following  form: 


dHn  V  Qn 

Ht~  “  "  TSSk)d  (Hn  -  “n-l1  +  TAtoT^ 


d¥ 

dt 


°ZPnz'nAxn  -  F 

%  _  Axn 


(5) 

(6) 


The  loop  subdivisions  are  given  in  Figure  6. 

EL  +  % 

For  section  1,  Q  «  for  section  5*  =«bV(L2_^ - hc) 

where  the  cooler  heat  transfer  coefficient  is  taken  to  be 
a  linear  function  of  velocity,  and  all  remaining  Q*s  are 
taken  equal  to  zero.  The  density  at  point  n  is  written  as 

Pn  -  *>«  +  %  <Kn  '  Hn  >  (7) 

o  o 

and  the  mean  density  of  the  nth  sub interval  is  set  equal 
to  the  average  of  the  densities  at  n  and  n-1. 

A  stable  solution  was  found  if  the  coefficient 
of  expansion  of  water  is  used  for  all  values  of  An 

oscillatory  solution  was  obtained.  Figure  7*  if  vaporization 


to  a  few  per  cent  qualify  were  assumed  in  the  top  section 
of  the  vertical  riser  (  =  -  9*500,  all  other  * 

-0.024  Tb/cu.ft . °F) . 
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Results 

Two  steady  equilibrium  modes  of  operation  were 
possible  when  the  pressure  at  one  point  in  the  loop  is 
held  constant  (surge  tank  E^  open  to  the  atmosphere). 

For  a  very  low  heat  input  the  water  temperature  in  the 
riser  never  exceeds  the  boiling  point  and  a  state  of 
stable  equilibrium  may  be  defined.  For  a  very  high  heat 
input,  the  entire  riser  contains  both  steam  and  water, 
and  a  maximum  flovj  rate  is  obtained.  Oscillatory  modes 
of  operation  result  for  the  intermediate  heat  input2.  An 
illustration  of  the  manner  in  which  the  period  and  amplitude 
depend  on  the  heat  input  is  given  in  Figures  8  to  12.  The 
period  and  amplitude  of  the  oscillations  are  determined  by 
the  mean  temperature  level  of  the  fluid  in  the  vertical 
riser. 

The  period  was  inversely  proportional  to  the 
mean  velocity  providing  some  steam  is  in  the  riser  at  all 
times.  When  the  system  does  not  contain  steam  during 
most  of  the  cycle,  the  period  is  considerably  longer  than 
that  predicted  by  the  extrapolation  of  the  higher  flow 
rate  periods. 

In  the  analogue  computer  problem,  =  3-3.6  Btu/sec, 
and  the  period  was  149  seconds,  which  is  about  21  seconds 
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less  than  the  expet imentally  observed  value  for  the  lowest 
flow  rate.  The  computed  period  would  be  expected  to  be 
less  than  the  observed  value,  since  the  density  function 
used  in  the  top  subinterval  in  the  riser  was  a  two-phase 
density  function.  Further,  the  density  function  used 
for  the  boiling  sub interval  was  linear  and  it  was  not 
possible  to  exclude  densities  greater  than  the  density  of 
saturated  water.  A  special  non-linear  element  is  required 
to  generate  a  density  function  with  the  correct  properties. 
As  a  result,  the  computed  mean  flow  rate  is  less  than  that 
which  could  actually  exist;  however,  the  flow  rate  curve 
and  enthalpy  curves  have  essentially  the  same  shape  as 
the  experimentally  determined  ones.  A  larger  computer 
is  required  for  solving  the  stability  problem  at  the  higher 
heat  fluxes. 


Conclusions 

A  natural  circulation  loop  can  be  made  unstable 
in  the  sense  that  a  small  displacement  from  the  equilibrium 
state  leads  to  undamped  oscillations.  Stable  operations 
result  when  the  fluid  temperature  in  the  riser  is  restricted 
to  values  less  than  the  boiling  point,  and  when  the  heat 
input  reaches  a  value  such  that  the  frictional  force 
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changes  wore  rapidly  than  the  driving  force.  The  theoreti¬ 
cal  treatment  of  an  open-ended  natural  circulation  system 
and  the  solution  of  a  simple  problem  using  an  analogue 
computer  have  lent  support  to  the  general  conclusions 
of  the  stability  analyses.  A  detailed  numerical  analyses 
is  in  progress  on  a  large  digital  computer. 
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NOTATION 

Equations  1  through  4 
A  =  cross-section  area  for  flow 

W  a=  total  mass  flow  rate 

x  8=  distance  along  streamline 

t  =  time 


P 

V 

H 

♦i 

k 

T 

P 


local  acceleration  of  gravity;  g  »  conversion  factor 


in  Newton's  law  of  motion 
=  slope  of  pipe  at  x 
=  wall  of  frictional  force 


fluid  density,  p  e  j  (AWPW  +  As  Ps) 

^  AwPw^w  +  VsVs  _  W 

volumetric  flow  rate.  V  S - ^ 

specific  fluid  enthalpy.  H  «  |pCAwPw\  +  AsPsHs^ 

VwVwHw  +  WsHs 


,Wn  *  VsVs 
wv 

=  thermal  conductivity 
=  fluid  temperature 
=  fluid  pressure 


<J>o  * 


m 


r 


=  radius,  rp  »  pipe  radius 
Q’  =  rate  at  which  heat  is  conducted  through  pipe. 
a(p)  =  an  empirically  defined  function  for  Eq.4  for  two-phase 
pressure  drop. 

Subscripts 

s  =  steam  phase;  w  =  liquid  water  phase 


For  the  computer  problem,  for  n  =  1,  2,  and  3,  H  =  l82Btu/lb, 

no 

Pn  -  59*781  lb/cu.ft . ,  #n  =  -  0.024  lb/cu.ft . °F;  n  =  4, 


=  182,  p  =*  37.5OO,  6fn  «  -9*500;  n  =  5  through  10, 
llo 


=  153,  Pn  -  60.477,  ^ 


-  0.024.  V0 
EL  +  Hj, 


=  3.16  Btu/sec ,  and  =  -  25.0V(~2~2“ 


0.0025  cu.ft./sec, 
-  89.6)  Btu/sec. 
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